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SHAPED SCINTILLATION DETECTOR SYSTEMS FOR 


MEASUREMENTS OF GAMMA RAY FLUX ANISOTROPY 

Abstract 

The detection efficiencies of cylindrical detectors for various gamma ray 
photon angular distributions have been studied in the energy range from .10 Mev 
to 15 Mev. These studies indicate that simple detector systems on small satellites 
can be used to measure flux anistropy of cosmic y-rays and the angular distri- 
bution of albedo gamma rays produced in planetary atmospheres. The results 
indicate that flat cylindrical detectors are most suitable for measuring flux 
anisotropy because of their angular response function. A general method for 
calculating detection efficiencies for such detectors is presented. 
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SHAPED SCINTILLATION DETECTOR SYSTEMS FOR 


MEASUREMENTS OF GAMMA RAY FLUX ANISOTROPY 
1. Introduction 

Recent observations of cosmic gamma rays and theoretical studies of their 
implications have indicated that observational data in the 0.3 MeV to 30 MeV 
energy range may provide important cosmological information about the early 
history of our universe. 1 " 7 Cosmic gamma-rays in this energy range are 
expected to be highly isotropic and thus of cosmological origin. A direct ob- 
servational test of the isotropy of the flux can be an extremely sensitive test of 
cosmological origin and such a measurement may provide a determination of 
the strength of galactic gamma-ray emission in 0.3 to 30 MeV energy region. 
Detailed spectral information in this energy range will provide important 
information as to the nature and origin of the gamma-radiation. 

It may be possible to study the structure of a planet’s atmosphere by 
measuring the energy and angular distribution of the albedo gamma-rays. Such 
measurements will also be considered in this paper. 

We present here the results of an investigation of the properties of cylin- 
drical detectors which can be used to optimize the design of simple systems for 
performing the following measurements: 

1. The gamma-ray flux anisotropy. 

2. The possible presence of point sources. 
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2. Definitions of Effective Cross Sectional Area and Intrinsic Detection 

Efficiency 

The response of a particular detector system can be evaluated in terms of 
a detection efficiency which depends on the geometry of the system and upon the 
spatial and angular distribution of the incident flux being measured. Because 
source to detector distances in astronomical measurements are very much 
greater than detector dimensions, spatial variations over the detector may be 
ignored. We will use a rectangular coordinate system (x 1 , y’, z') fixed in the 
detector to describe positions within the volume of the detector system. We 
will denote a spatially fixed coordinate system by the unprimed rectangular 
coordinates (x, y, z) and use (b, 6) for the latitude and longitude angles. When 
one considers the application of the results derived in this section to gamma 
ray astronomy, the most convenient unprimed system will be the galactic system. 
The angles (b, -i ) will be used to describe the direction of the incident photons. 
The orientation of the detector system will be specified by giving the latitude 
and longitude of the z' and x 1 axes denoted by b^, , £ , and b^, , •£ , , respectively. 

The geometry of a general situation is given in Figure 1. The x' and y' axes 
are not included in order to show more clearly the various angles and a photon 
flux from direction (b, <t) entering the detector at point A (x' f y' f z’). The 
rectilinear path through the crystal along the incident photon direction is denoted 
by AB with a distance £ . We will denote the photon flux which is differential in 
angle and energy by the quantity 
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(1) 


f(E, b, t) = ^ dE 

di2 

where 

dH = d(sin b) dl. (2) 

We can then compute C, the total number of photons per unit time and energy 
which will interact with the detector, by 


C (E, b z , , '^ z ' > b x < > » ) 


= 1 f 

Jfl Js' 


f (E, b, l) 


( 3 ) 


{1 - exp [-p.(E) £ (b z ,, -i z ,, b x , x\ y f , z', b, 1} } dA^ dQ, 


where dA ± is the element of detector surface area perpendicular to the direction 
(b, -C) at the point of entry A (x’, y’, z*) and S’ represents the Integration over 
the surface of the detector where photons enter. The quantity /j. is the linear 
absorption coefficient which is energy dependent and represents the sum of the 
photoelectric, Compton, and pair production interactions. The quantity within 
the braces represents the probability that the photon will have at least one inter- 
action within the crystal. Multiple interactions occur in times shorter than the 
decay of the scintillation produced in the detector and thus appear as a single 
interaction. These multiple interactions will change the amount of energy 
deposited in the crystal but not intrinsic efficiency. 

The effective detector cross sectional area is defined as: 

"-ftr 

Ja 
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and the intrinsic detection efficiency as: 


f f fdDdA ± ^ 

Jn J s ' 

The variability of e as a function of energy for three special flux distribution 
is illustrated in Figure 2 for a relatively symmetric crystal 10.16 cm long x 
10.16 cm in diameter. More extreme shapes most likely will show greater 
differences because of increased proportions of short path lengths for an isotropic 
flux. Therefore, the quantities t? and e are not useful in the case of an unknown 
photon distribution. They are also difficult to calculate for irregular detector 
geometries, but special cylindrical cases which are useful will be examined. 

3. Response of a Cylindrical Detector to a Parallel Beam 

We now specialize the geometry of the crystal to a right cylinder of height, 

H, and radius, R, with z f as the axis of symmetry and consider the flux distri- 
bution to be a plane wave coming from direction (b 0 , ^ 0 ), We do the plane wave 
case first since the results for a planar and an isotropic flux can be calculated 
in terms of the plane wave case. For a plane wave 

f = P(E) 5 (sin b - sin b„") S (•£ - 'Ll (61 

v ' \j ' ' ' 

Since we have cylindrical symmetry, the only angle of importance is the one 
between z' and the direction of the photon beam. We will denote this angle by 9 
and it is given in terms of (b 0 , ) and (b , , , ) by 

cos 8 = cos b 0 cos b z , cos (Z Q - ^ z ,) + sin b Q sin b £ , , (7) 
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and P(E) has units of photons/cm 2 -sec-MeV. In this case we obtain 


V(E, #) 


f r 77 / 2 



-2RHSin#-«!l - ( cos/3 

exp-f 2 ^ Rcos ^ 

d/3 

l Jo 

\ s i n # /_ 

J 



for tan 8 2 2R/H, and 


V (E, 8) = 2RH sin# 



cos/3 


exp 


f 2^R co s /?\ 
\ s i n 6 ) 


d/3 + sina 



+ 7T R 2 


cos 




p7T/2 

cos/3 



2,aR cos 
s i n 6 





4R . Q n 
s 1 n 6 co s 9 




d/3 -sina 



for tan# <: 2R/H, where Cos a ~ H/2R tan # and cos fi = £/2R. 
The intrinsic detection efficiency is then given by 


s 0 ) 


y*) 

2RH sin # + 77R 2 cos # 


( 10 ) 
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This case is useful when one is dealing with a single distant point source 
since the detector counting rate will be given by 

C ps (E, 6?) = V E ' *> P (E) <U) 

for a given energy photon E and by 

CD 

C PS ^ E t’^ = J % p ( E > dE (12) 

for all photons above a threshold energy E T . 

4. Response of a Cylindrical Detector to a Planar Flux 

For simplicity of calculation we consider the planar or disc like flux to lie 
in the plane b = 0. In this case 

f = J (E) 8 (s in b) (13) 

and J has units of photons/cm? -sec-MeV-rad. We wish to consider only two 
orientations of the cylindrical detector. Case 1: The z' axis lies in the plane 
b = 0 (i.e., z * z’ = 0). Case 2: The z' axis is perpendicular to this plane (i.e., 
z • z f = | z | |z |). For Case 1, we find from (7) that 

6 = •i — K> > 

z 

and we can choose l z >= 0 without loss of generality. Then we have, by performing 
the integration in (3) , 

r 2n 

Cj (E, b 2 , = 0) = J(E) 1? P (E, 9) do (14) 
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Since 


and 


we obtain 


and 


J 


f dQ = 2ttJ (E) 



fdndA ± = 4[2RH + n R 2 ] J(E) 


7? (E, b 2 i - 0) = - 

J 77 


«j(E, b t , = 0) = 



7 p (E, d) dff 


77 p (E, 6) dd 


[2RH + ttR2] 


(15) 


(16) 


(17) 


(18) 


For certain purposes, it is useful to define an effective geometric factor in this 
case by 


Gj(E, b 2 , = 0) = 
2 

We note that Gj has units of cm -rad. 
6 - 77 /2 and it follows that 


J (E) 


r>n/ 2 

= 4 J 

•^0 


7J p (E, d) dd 


(19) 


For the 2nd case we have from (7) that 


b 2 , - 77/2) = 2 ttJ(E) 77 p (E, 77/2), 

(20) 

(E, b g , = 77/2) = V p (E, 77/2), 

(21) 

fdndA ± = 4 t7RHJ(E), 

JnJ S ' 

(22) 
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(23) 


€j(E, b z » -tt/2) 


V E > tt/2) 
2RH 


and 


G J (E, b z / = n/2) = 2 ttt] p (E, tt/2). 


(24) 


5. Response of a Cylindrical Detector to an Isotropic Flux 
In the case of an isotropic flux 

f (b, -C , E) = 1(E), 


(25) 


where 1(E) has units of photons/cm2-sec-MeV-ster. For simplicity of calculation 
and without any loss of generality we take the (x 1 , y 1 , z') system exactly coincident 
with the (x, y, z) system. Then we find from (7) that 8 = tt /2-b 0 and 


»tt/2 

I ^ 

^0 


Cj(E) = 4 ttI(E) 7] p (E, 0) d(c o s 0 ), 


7?! (E) = 1 


V p (E, 8) d(cos 8 ), 


(26) 


(27) 



f df)dA A = 2 tt 2 RI(H + R), 


(28) 


and 


e x (E) 


2 

tt-R(H + R) 



7j p (E, 8) d(cos 8). 


(29) 
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In this case we also define an effective geometrical factor by 


C (E) 

G I (E) = = 47T J r) p (E, 0) d C cos d) (30) 

and G r has units of cm 2 -ster. 

6. Measurement of Cosmic Gamma Ray Flux Anisotropy with a Cylindrical 

Crystal 

The effective geometry factors, G y (E, O), Gj (E, tt/2) and G r (E) can be 
used to show how the anisotropy of the cosmic gamma ray flux in the energy 
range 1-10 MeV can be measured or have a sensitive upper limit placed on it. 

If the flux is not isotropic, the most likely source of excess photons will be 
from the galactic plane. We consider a spinning satellite placed in an orbit such 
that it spends a considerable fraction of its time away from the particle radiation 
environment near the earth and has its spin axis in the galactic plane. Spallation 
effects produced by cosmic rays in Nal and, C si crystals are not considered, 
although they have been shown (Reference 11, 12) to be an important background 
in determining the true cosmic gamma ray flux. However, these will not interfere 
with the measurement of the flux from the galactic plane by this method, but will 
increase the background and the error in the determination of anisotropy. 

A flat cylindrical (disc) crystal placed in the spinning satellite with a diameter 
lying in the spin axis will result in a motion such that the z’ axis will spend some 
of its time perpendicular to the galactic plane and some of its time in the galactic 
plane. We show later that a disc is better than a rod for a given mass detector. 
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We denote the total amount of time spent in each position to be T; then the total 
counts per unit energy will be for the case of z* perpendicular to the galactic 
plane 


qcE) = 


i 


1(E) Gj(E) + J (E) G 



and for the case z f parallel to the galactic plane 


(31) 


C„ (E) = [I (E) Gj (E) + J (E) Gj (E, 0)] ■ T- (32) 

These can be solved for J(E) and 1(E) to obtain 


J = 


and if we assume that I > > J 


(C|| c x > 

O)- 0 ; 


H] 


(33) 


I 



(C„ -C^GjCE, 0) 1 C|| ^ 

GjCE, 0) - Gj ^E, |)| G i T G i T 


(34) 


We know from experimental evidence that J is small in the energy range of 
interest and therefore C ± = . We will define a geometrical sensitivity factor 

for a monoenergetic flux by the quantity 


S' = 


G ; (E, 0) - G 



Gj (E) 


(35) 


and the effective anisotropy by the quantity 

s 1 

4ttI 21 


( 36 ) 
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Since I - W ^ Background, we cannot measure true anisotropy 
method. The effective anisotropy will be less than the true amsotropy. Wi 

this understanding, we can write 


c ll ' c i ^ 

S = ^q~ S' 


( 37 ) 


The statistical error 


(i.e. standard deviation) in 8 is 


= 2S 5 


'(C 11 + q ) 1/2 (q-c^q 


1/2 


q 


.L&. 

S' V 2C l 


(38) 


. 1 . • +Vinf fh<=* minimum m6EisuF6£ibl® 

. „„ r . c Then if we use the criterion that the munm 
since C n - ‘-'i • . .. . 

qu antity must he 3 times the standard a—, the — value of 
can be observed is 


.3_£T. 

s' '2q 

of I and T, we use (34) to obtain 


(39) 


S min = “ S r »'2q 


Now to write this in terms 

8 = a /XT = l^ 

°mir\ g' y2IG I T S 2 IT 


(40) 


where 


s = s' /q = 


3j(E, 0) -0, (e.-|) 




(41) 


Thus for a fixed value of I, one can decrease 
Of 8 by making | Si as large as possible and by increasing T. 


the minimum measure able value 
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In the case of an integral spectrum we can follow a similar analysis by 
starting with (31) and (32) and integrating over E. 

C l( >E)=[7 * (E) Gj (E) dE + | J(E)G J ^E,^dE T ( 42 > 

C |( (> E) = jj£ 1(E) Gj(E) dE + J J(E) Gj(E, 0) dE T ( 43 ) 

Another useful quantity can be obtained from rewriting equations (42) and 
(43) in the following way 


C i (>E) = Gj I( E) dE + Gj f - ) T (E) dE 


q i (>E) = G I 1(E) dE + G 


i (0> f 


J(E) dE, 


where 


T 

•'E 


1(E) dE 


1(E) dE 


J(E) dE 


77 \ ~E 


J(E) dE 
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and 


G J (0) = 


r 


Gj (E, 0) J(E) dE 


(46) 


J 


J(E) dE 


These average values of the G's can then be used to determine the sensitivity 
factor S and S' for particular spectral distributions of interest. 

To illustrate some points we consider the case of an opaque detector i.e., 

M - 00 . In this case, we see from equation (8) and (9) that 

r] p (8) = 2RH Sin 9 + ttR 2 Cos 6. 

This means that 


Gj (0) = 4R [2H + 77 R] , 


(47) 




= 477 RH, 


(48) 


G r = 2 t7 2 R(H + R), 


(49) 


and 


2/2 R 

j-l 

^77 -2\ 1 
n 

\ 77 / _ 

/ a + 1 


(50) 


where a - H/R. 
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If S = 0, the detector will have no sensitivity to flux anisotropy for the two 
orientations considered above. From equation (50), S will be zero for 

a = 2.75 

TT - 2 

or 

H - 2.7 5 R 

Thus for the limiting case of an opaque crystal, when the height becomes some- 
what greater than the crystal diameter, the ability to detect flux anisotropy by a 
rotating crystal becomes very small. 

The sensitivity factor S for Csl is shown as a function of height for a number 
of crystal radii in Figure 3 for ^ - 4 (^ .15 MeV), in Figure 4 for ^ - .4 (~.56 
MeV), and in Figure 5 for ^ = .165 (- 3 MeV to 7 MeV). 

The choice of crystal size for maximum sensitivity to the proposed flux 
anisotropy can be determined from such sensitivity curves. The curves tend to 
indicate maxima for the higher energy gamma radiation (Figure 4 and 5). 
Although not shown in the figures, it has been found that for a given volume 
(mass) of crystal, it is better in terms of sensitivity to use a disc (2R > H) 
rather than rod (H > 2R). 

As an example, the results of calculations of the minimum measure able 8 
for a 2 inch thick by 8 inch diameter Csl crystal are presented. For this case 
the isotropic gamma ray spectrum (approximated by I = .024E” 1 - 7 ) between .3 
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and 10 MeV as measured by one experiment on Apollo 15 (Reference 7) was used. 
A power law spectrum with a variable exponent was used for the planar flux. 

The results are not sensitive to this exponent for spectra harder than an index 
of minus 2. G r , (O), and Gj (77/ 2) were calculated for these spectral distri- 
butions. A flight time of 90 days was assumed with a 10 percent duty time in 
each orientation with respect to the planar flux. Using equations (40) and (44) 
yields the minimum effective anisotropy which can be measured by the system 
of .00038. Using equation (36) and the isotropic flux yield a total planar intensity 
in the .3 to 10 MeV range of 5.5 x 10 _s photons/cm 2 -rad -sec for differential 
energy spectrum with a minus 1 exponent. 

7. Measurement of Sources in the Galactic Plane 

One can examine the use of this simple system for the detection of sources 
in the galactic plane. Essentially, the detection efficiency, rj p , (E, 8 ), in equa- 
tions 8 and 9 reflects the detector's response as a function of phase angle 6 
(i.e., the angle between the direction of the source and the cylindrical crystal 
axis). Figure 6 is a polar plot of the ratio 


X (E. £) 



as a function of phase angle 6 for a 40.64 cm x 2.54 cm Csl right cylindrical 
crystal. The calculation was made for approximately 4 MeV gamma rays. 
Maximum sensitivity is obtained with the face of the crystal pointing toward the 
source {0 = 0). As one rotates the crystal the count rate decreases until 6 -n/2. 
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The count rate decreases by a factor of 4.5, with the greatest change oceuring 
near minimum. Thus, a source would be better noticed by the minimum rather 
than the maximum in count rate. We can define an aperture or angular resolution 
as the width in degrees at half the counting rate. The angular resolution of this 
system would be ~ 18°. 

The sensitivity of the system to the detection of point sources is obtained by 
calculating the expected count for 8 - 0 and 0 = 77/2 with respect to the source 
in the presence of the isotropic flux considered in the previous section. The 
total expected counts, C i , for 8-0 would be 


C = T f l7 ?I dE+ f 

Je Je 


P77 (E, 0) dE 


where P is the parallel beam flux due to the source in photons/cm 2 -sec-MeV. 
The total expected counts for 9 = 77-/2 would be 


c 2 = T 


J I^dEt [p^,(e, |)dE 


We now define 



I 

I 

J, 


IdE, 

PdE, 
Iijj dE 


I 
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and 


%( 9 ) = 


s. 


Pv p (0) dE 


Then the point source strength averaged over energy can be found from the 
difference in counting rate and known average crystal properties by 


C - C 

'-'i <~ 2 




The standard deviation, <7 of P for weak sources is 

p 


(51) 


- / 2 I 17 t • 


(52) 


For detectability, again we require that 


Assuming an integrated omnidirectional flux, T of 7.2 x 10“ 2 photons/cm 2 -sec- 
ster and a measurement time of T = 2.6 x 10 s sec in each position (similar to 
the situation discussed previously), P min > 4 x 10~ 4 photons/cm 2 -sec, could 
theoretically be detected. Galactic sources such as the Crab Nebula are about 
an order of magnitude greater than this. Systematic errors can, of course, 
decrease this sensitivity. 
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8. Measurement o f Atmospheric Albedo Gamma Hay Distribution 

This methodology of reducing calculations to integrations over the parallel 
beam case can be applied to more general flux distributions than planar or point 
sources. An example is the atmospheric albedo gamma ray flux. For this case 
the space coordinates are such that the z axis extends from the detector to the 
nadir. The flux is considered to be independent of •i. Calculations are done for 
a monoenergetic flux. 

The albedo is much greater than the cosmic flux, so the isotropic background 
is not important in this case. The main concern here is the detailed variation of 
the counting rate with detector orientation. Equations (8) and (9) are again used 
to define y . Using equation (7) to define & in terms of the flux and detector 
orientation angles, the counting rate for a given orientation can be written as. 


C(E, b z 0 


r 2n r max 
Jtt 



b) K(b) dbd{X - % z <)- 


The independence of C to l z <, is a result of the flux symmetry. The relative 
counting rates for a specific size crystal are given in Figure 7 for various 
assumed angular distributions. Two different limb angles corresponding to 
different altitudes for the detector are shown. At higher energies Stecker 


(Reference 8) has proposed a model for and calculated the scanning angle distri- 
bution of albedo gamma rays. This distribution was found by Stecker to be well 


approximated by a secant up to about 50°. It was found that a disc shaped 


detector gave a greater variation in counting rate for this model than an equal 
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volume (mass) rod shaped detector. In addition the disc gives a better differen- 
tiation between different exponents of the secant. It can be seen that if the limb 
angle is known the departure from a secant distribution at the medium energies 
(3-10 MeV) can be determined with a disc shaped detector. 

9. Summary 

The cases which have been described involve either determining the intensity 
of a particular flux in a predetermined background or des criminating among a 
limited number of flux distributions. The general formulae can be used to de- 
scribe any flux distributions which can be reduced to integrations over parallel 
beams if an iterative scheme is used. Iteration is necessary because the quan- 
tities V and € used to calculate the flux from the counting rate are themselves 
dependent on the flux distribution. After the first analysis of the data the effi- 
ciencies can be redetermined in order to calculate the flux distribution again, 
and so forth. 

We feel, therefore, that a flat cylindrical crystal in a small satellite to keep 
local background problems to a minimum is a viable system for studying gamma 
ray flux distributions from orbit. 
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Figure 2- intrinsic Efficiencies as a Function of Energy for 10.16 cm x 10.16 cm (4 x 4 ) Cylindrical Crystals. 








SENSITIVITY = (G 2 ~ G x )/ VG^ 



figure 5. The Sensitivity Factor S for Csl as a function of height for a number of crystal radii for /x -.165 (~ 3 MeV). 



Np(E,S) 

Figure 6. Polar plot of the relative effective detector cross sectional area ^ ^ 77/2) as a 
function of 6 for a 40.£2 cm. x 2.54 cm. thick cylindrical Csl crystal and 4 MeV protons. 
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Figure 7* Calculated Counting rate variation relative to maximum as a function of detector orientation in atmospheric 

albedo; 8.89 cm rad x 2.5 cm thick crystal, 1.5 MeV gamma rays. 




